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5-H ' Abstract. In this paper we derive structure theorems that characterize the 

tJn' spaces of linear and non-linear differential operators that preserve finite di- 

■^^ , mensional subspaces generated by polynomials in one or several variables. By 

1 means of the useful concept of deficiency, we can write explicit basis for these 

>D ' spaces of differential operators. In the case of linear operators, these results 

Cn ' apply to the theory of quasi-exact solvability in quantum mechanics, specially 

in the multivariate case where the Lie algebraic approach is harder to apply. 
In the case of non-linear operators, the structure theorems in this paper can 
r^) ' be applied to the method of finding special solutions of non-linear evolution 

• ' equations by nonlinear separation of variables. 
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^ . 1. Introduction 

O 

C^^ ' It is a fact that the Schrodinger operators whose point spectrum, or at least 

^P , part of it, can be computed algebraically are often related to differential operators 

admitting invariant spaces of polynomials. Lie algebras have played a unifying role 
^^ , in this area, because many of these polynomial spaces turn out to be irreducible 

f^ ' modules for a faithful representation of a finite-dimensional Lie algebra by means of 

first-order differential operators. The classical theory of quasi-exactly solvable po- 
tentials has thus been built on the assumption that the exactly solvable Schrodinger 
r^ , operator under study should be expressible as a quadratic element in the universal 

enveloping algebra of a finite-dimensional Lie algebra of first-order differential op- 
erators, admitting an explicitly computable invariant subspace of square-integrable 
functions, or a complete infinite flag thereof [1-4]. Burnside's Theorem serves as 



> 

X, 

JH , a strong argument in favor the Lie algebraic approach since it implies that any 

endomorphism of an irreducible module for a Lie algebra can be represented as 
a polynomial in the generators of the algebra. However, recent results show that 
the Lie algebraic approach suffers from various limitations that reduce severely its 
applicability: 

(1) In the case of polynomial subspaces in one variable, the Lie algebraic ap- 
proach can only be applied to find the differential operators that leave the 
polynomial space Vn = (1, z, z^, . . . , z") invariant, but it cannot character- 
ize the set of differential operators that map Vn into Vm C Vn with m < n. 
This simple problem has motivated the important notion of deficiency used 
throughout this paper, and applications of it can be found in the construc- 
tion of solvable classical many-body problems by considering the motion of 
zeros of polynomials whose coefficients evolve in a controlled manner [5] . 

1 
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(2) Other subspaces generated by polynomials exist which are not isomorphic 
to Vn- The Lie algebraic approach cannot be applied in these cases, as was 
already noted by Post and Turbiner, who characterized the spaces of linear 
differential operators which preserve polynomial subspaces in one variable 
generated by monomials. In their work [6] they solved this problem with no 
reference to Lie algebras. The case of a general space spanned by polynomi- 
als — referred to as the generalized Bochner problem in [6] — remains still 
open. Somewhat surprisingly this direct approach has not been pursued 
until very recently, where it has been shown that the class of quasi-exactly 
solvable potentials is larger than the Lie-algebraic class [7] , and that even in 
Lie- algebraic potentials other non-s[2 algcbraizations exist which allow to 
obtain more levels from the energy spectrum of the Hamiltonian, [8]. The 
existence of differential operators that preserve a general polynomial space 
and cannot be expressed as quadratic combinations of the generators of 5(2 
is not in contradiction with Burnside's theorem, since a general polynomial 
space is not the carrier space for an irreducible representation of SI2. 

(3) In the case of multi-variable polynomial subspaces, the problem of char- 
acterizing the set of linear differential operators that leave these spaces 
invariant becomes untractable in the Lie algebraic approach. The reason is 
that, contrary to the single variable case where essentially s[2 is the only 
algebra of first order differential operators with finite dimensional repre- 
sentations, in more variables many more algebras exist. But moreover, the 
characterization of second order operators as quadratic combinations of the 
generators of these algebras requires an extensive analysis of the syzygies 
corresponding to the primitive ideals associated to the irreducible represen- 
tations [9] . These problems are entirely bypassed in the direct approach, as 
shown in Section 3 of this paper, where a simple characterization is given 
for the set of linear differential operators of any given order r that leave 
the simplicial module 

Vn — spanjx'j^ . . . x^ \ii + ■ ■ ■ + in <n} 

invariant. Our results coincide with the formulas for the special case N — 2 
and r ^ 2 derived in [9] using the Lie algebraic method. 

It has now become clear that the connection to Lie algebras is not an essential fea- 
ture of exact or partial solvability. Our goal in the first two Sections of this paper is 
to present a direct method to characterize linear differential operators with invari- 
ant polynomial subspaces which is simpler and more powerful than the Lie algebraic 
approach. We restrict in this paper to the simplest case of polynomial subspaces, 
namely the simplicial modules Vn, the case of general polynomial subspaces shall 
be treated elsewhere. For any number of variables A'^, we provide an explicit ba- 
sis for the space of linear differential operators of any order r that map Pn into 
Vm C Pn with m < n. It should be stressed that although these results allow to 
construct many differential operators with invariant finite dimensional polynomial 
subspaces, in general it is not known whether a transformations exists that puts 
the operator in Schrodinger form. Therefore the results in Sections 2 and 3 are only 
a first necessary step in the theory of higher dimensional quasi-exact solvability. A 
general theory would need to face the difficulties of the equivalence problem [10]. 
Despite this fact, it is worth mentioning that a few examples of partially solvable 
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multi-dimensional Hamiltonians exist [11-13], which are mostly extensions of the 
Calogero-Sutherland class. 

This paper also addresses the study of nonlinear differential operators with poly- 
nomial nonlinearities which possess invariant polynomial subspaces. The motiva- 
tion for this study is twofold: 

(1) In [14], the important concept of operator duals is introduced. Given a 
finite dimensional space of functions J-'n = span{/i, . . . , /„} which are re- 
quired to satisfy certain regularity conditions, the operator duals are linear 
differential operators defined by the relations 

These operators are used in the reduction of non-linear evolution equations 
to dynamical systems by a method of non-linear separation of variables. 
In [14] the existence of these differential operators is proved together with 
results on the regularity of the coefficients. The proof is constructive and 
therefore given any space !Fn whose basis elements satisfy the required 
regularity conditions, non-linear evolution equations can be written which 
have solutions in the space Tn, i-e. special solutions exist of the form 
u{t,x) = X]"=i Ci(i)/i(2;) where the coefficients Ci{t) satisfy a system of 
coupled non-linear ODEs. However, the order of the operator duals is 
precisely the dimension of the space JF„. The price to pay in order to have 
invariant spaces of high dimension is that the order of the resulting evolution 
equations grows with the dimension of the space. For dimensions n larger 
than six this reduces the applicative interest of the resulting equations. 
The motivation coming from [14] is to construct operators that generalize 
the operator duals, so that the order of the resulting equation and the 
dimension of the invariant subspace can be independently chosen. In the 
case of polynomial subspaces P„ , the generalization of the operator duals 
to arbitrary order r < n are the so called maximal deficiency operators 
introduced in Section 2. 

(2) The second motivation comes also from the results on non-linear separa- 
tion of variables by King, Galaktionov and Svirshchevskii [15-19]. In these 
papers special interest is given to translation-invariant evolution equations 
with quadratic non-linearities which admit solutions via non-linear separa- 
tion of variables. From a physical context, applications are found in nonlin- 
ear diffusion and thin film equations. In this paper we extend these results 
by providing a comprehensive structure theory for autonomous nonlinear 
operators that preserve a polynomial space Vn- 

Our paper is organized as follows. In Section 2, we present our direct approach 
in the case of linear differential operators in one variable. Besides the order of a 
differential operator, we introduce two key invariants which can be freely specified, 
which are the degree and the deficiency of the operator relative to a polynomial 
space Vn- The order, degree and deficiency are shown to specify the operator 
uniquely up to scaling by a constant. The operators of maximal deficiency generalize 
the operator duals of [14] to any order lower than n. We give an explicit basis for the 
space of operators of given order and deficiency. Section 3 is concerned with linear 
differential operators in N variables, where all the results of Section 2 are shown to 
extend to the case of simplicial modules, that is multivariate polynomials of total 
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degree bounded by a given integer. Section 4 is concerned with non-linear operators 
preserving polynomial modules, where we give an explicit decomposition theorem 
for the most general non-linear operator with polynomial coefficients preserving a 
siniplicial module. Section 5 studies the deficiency concept for non-linear operators 
with polynomial nonlinearities. Section 6 concentrates on non- linear operators that 
are translation invariant (autonomous) deriving also structure theorems for this 
class. On Section 7 the application of these results to non-linear separation of 
variables is discussed while some explicit formulas are given for quadratically non- 
linear autonomous operators in Appendix A. 

2. Linear operators in one variable 

In this section, we consider the class of scalar linear differential operators on 
the real line, with polynomial coefficients. We are interested in the subclass of 
operators which have a definite order, degree and deficiency. These quantities are 
defined as follows. The order of 

(1) L = j2a^i^)D\ D:=^, 

i=0 

is as usual the largest r for which the coefficient ar{x) is not identically zero. We 
say that L is of degree d G Z if for all j E N, there exists Cj ^ E M. not all zero 
such that 

(2) L[x^]^cjx^+'^. 

In order to define the deficiency, we fix n G N and consider the vector space 

(3) Pn^spa,n{l,x,...,x"^} 

of polynomials in x of degree less than or equal to n. We say that L has deficiency 
m G Z relative to Vn if 

(4) LVn C Vn-m, but LVn t Vn^m-l- 

Let Cr^m denote the set of linear differential operators with polynomial coefficients, 
of order less than or equal to r and of deficiency greater than or equal to m relative 
to Vn- Again, we emphasize that the notion of operator deficiency only makes 
sense relative to a particular n. Most of our discussion will be carried out with the 
assumption that the n in Vn has been fixed. As such, we will often omit the n in 
our terminology and notation, simply speak of the deficiency of an operator, and 
write Cr,m instead of Cr\m- 

Proposition 1. The set Cr^m is a subspace of the vector space of all linear differ- 
ential operators, i.e., it is closed under linear combinations. 

Proof Given linear operators L, L' G Cr.nn the order of any linear combination of 
L and L' is less than or equal to r. Similarly, the deficiency of a linear combination 
is greater than or equal to m. D 

Proposition 2. The deficiency of a non-zero linear operator cannot exceed its 
order. 

Proof Suppose that L G Cr.m is a linear operator such that m > r. The operator 
jjn-m+ij^ annihilates Pn, but has order less than n + 1. This is impossible. D 
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Trivial examples of operators with given order, degree and deficiency are given 
by the operator £)*, which has order i, degree —i and deficiency i, and the multipli- 
cation operator x^ has order zero, degree j and deficiency — j. These operators do 
not depend on the degree n of the polynomial space Vn- A more significant exam- 
ple, which depends explicitly on n, r, to, d with < m < r < n and —m < d < r~m 
is the operator 

(5) Lrmd '■= x^ (n — j — xD )kD\ i = r — m, j = r — d — m, k = d + m 
where we have introduced the Pochhammer operator 

(6) {n-xD)k:^ {-!)''{ xD ~ n){xD - {n - 1)) ■ ■ ■ {xD - {n - k + 1)). 
A basic result is the following: 

Proposition 3. The operator Lrmd has order r, degree d and deficiency m. 

Proof Since the fc-fold composition appearing in the right-hand-side of © anni- 
hilates the monomials x", . . . ,a;"~'^+^, the operator in (jSJ has order k, degree zero 
and deficiency k. It follows that the operator {n — j — xD)kD^ has order r = k + j, 
deficiency m ^ k + j and degree d ~ —j. By left-multiplying by a monomial of x 
we raise the degree and lower the deficiency, so that Lrmd has order r, deficiency 
TO, and degree d. D 

Following Proposition 12 it is helpful to refer to an operator whose deficiency 
equals its order, as a maximal deficiency operator. To this end, we will use the 
symbol 



^L -^- 



(7) Krj ^ -, —Lr^r.-j = j —{n-'j-xD)r-jD\ < j < r < n, 



to denote the operators of maximal deficiency. The normalization constant of >^_ ■>, 
will be useful in later formulas. These operators have a number of interesting 
properties, and play a key role in our theory. 

Proposition 4. Up to multiplication by a non-zero real constant, the operator Krj 
is the unique r**^ order maximal deficiency operator with polynomial coefficients 
having degree d — —j , where < j < r. 

Proof Let L be an r^^ order maximal deficiency operator with polynomial coeffi- 
cients and having degree d. Since L has polynomial coefficients, d > —r. As well, 
L maps x'' to a multiple of x''^'^. Hence, 

L[x''] =0, n-r-d<k<n. 

However, a non-zero r^'^-order operator can annihilate at most r monomials. Thus, 
we have established that < j < r, where j = —d. 

The leading order term of both L and Krj is a multiple of x'^^^D^ . Hence, 
there exists an a £ M such that the order of L — aKrj is less than r. However, the 
deficiency of L — aKrj is greater than, or equal to r. Therefore, by Proposition |21 
we have L = aKrj. D 

Proposition 5. For every < j < r < n, we have 

(8) (r - j) ! Krj = {n - i - xD)r-j D' ^ D\n - xD )r_j . 
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Proof The operators {n — j — xD )r~j D^ and D\n~ xD )r-j both have order r, 
deficiency r and degree —j. By Proposition 01 they differ by a scalar multiple. By 
comparing the coefficients of the leading order, we see that these two operators are 
actually equal. D 

Proposition 6. For a fixed r, the operators K^j are recursively defined by 

(9) [D,Krj]^-Kr,, + l, Krr = D'-. 

Proof Setting k = r — j, we have, by Proposition[3 

[D, Kr,] = -^— {D^+\n - xD)k - [n - j - xD)kD^+^) 

^ D^+\{n - xD)k - {n + \ - xD)k) 



{r - j) 



^ D^+\n~xD)k-i 



(10) ^n=Y.i-^r('':l:')i^'^' 



= -Kr^j + l- □ 

We can also expand the maximal deficiency operators into an operator sum. 
Proposition 7. We have 

, .. . , IT 

fc=0 ^ ' 

Proof Let K^j denote the right hand side of pU|) . A direct calculation shows that 

Hence, by Proposition|Hl Krj and Krj have the same recursive definition. Therefore, 
the two operators are equal. D 

Proposition 8. We have dimCr^r = r + 1. Indeed, the operators Krj, j = 
0, 1, . . . , r form a basis for the vector space Cr^r- 

Proof Let L e Cr^r be given. Decomposing L into operator monomials of homo- 
geneous degree, we have by Proposition [5] and the fact that Vn is generated by 
monomials, 

r 

L ^ 2_^aj Krj ■ 
i=o 

Since the Krj are linearly independent, we conclude that dim Cr^r — r + 1. D 

The following corollary follows immediately. 

Proposition 9. We have 

(11) Cr.m — 'Pr-m(S> Cr.r, with dim £,..„!:= (r + 1) (r — TO + 1). 

In particular, the set of operators of order r or less that map Vn to itself is of 
dimension given by [r + 1)^ = dimgl(r + 1,M). 

We conclude this section with a simple example: 
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Example 1. A basis for the vector space £2,2 of operators of order two and 
deficiency two is given by 

K22 = D^, 

K21 = {xD-n + l)D = xD'^ + {l-n)D, 

K2a = {xD-n){xD-n+l) = x'^D'^ + 2{l-n)xD + n{n-l). 

3. Linear operators in several variables 

The results of the preceding section extend readily to the case of linear operators 
in A'^ variables (xi, . . . , xn). We consider the vector space 

(12) T'^ = {x\' ...x'^ \n + --- + iN<n}, 

of polynomials of degree n in iV variables, of dimension ( ). We shall use the 
standard multi-index notation whereby given a multi-index / = (ii . . .ztv) G N^, 
we let 

(13) x^:^x\^...x^, Dj:^^---^. 

ox^ ox^ 

The notions of order, degree and deficiency are defined similarly to the single- 
variable case. Thus, an operator 

r 

(14) -^ = X! "-f (^1' • ■ • ' xn)Di, 

|/|=o 

will be of degree d £ Z if for almost all / e N^, there exists cjj ^ E (M^)^ such 
that 

(15) L[x']= J2 c/ja;'^. 

\J\ = \I\+d 

In order to define the deficiency, we fix again n G N, and say that L has deficiency 
771 e Z relative to V^ if 

(16) LV^ C T'^.™, but LV^ t r^_^_,. 

As in the single variable case, it is easy to see that the deficiency of an operator can- 
not exceed its order. For example, the operator Dj has order |/|, degree — 1/|, and 
deficiency |/|, and the operator x^ has order zero, degree |/|, and deficiency — 1/|. 
Again, a more significant example is obtained by introducing the Euler operator 

(17) ^-E-.^j. 

the Pochhammer operator 

(18) (n - E)k -.^ (-1)''{E - n){E - {n - 1)) . . . {E - (n - k + 1)), 
and considering the operator 

(19) x'in~\J\-E)kDj. 

This operator has order | J| + A:, degree |/| — \J\ and deficiency | J| + A: — |/|. Again, 
if we let Cr,m denote the vector space of linear differential operators in N variables 
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(xi, . . . ,xn)j with polynomial coefBcients, of order r and deficiency m, then we 
have 

20 dim Cr,„, = { 

\ r — m J \ r 

or equivalcntly 

This formula is in agreement with the result obtained in [9] for operators in two 
variables, in which it was proved that 

(22) dim Crfi = 



2+ ^' 



r 

The proof given in [9] was less direct and required an analysis of the syzygies defined 
by the primitive ideals associated to the irreducible representations of sl(3,M). 

It is easy to see that in contrast with the single variable case, the order, degree 
and deficiency are not sufficient to characterize an operator uniquely up to a non- 
zero factor. We have: 

Proposition 10. A basis for the vector space of linear differential operators with 
polynomial coefficients, of order r, degree d and deficiency m, in N variables is 
given by 

(23) x\n - \J\ - E)kD J , \I\ = r - m, \J\= r - d - m, k = d + m. 

This vector space is thus of dimension 

,„,, fN + r-m-l\fN + r-m-d-l'^ 

\ r — m J \ r — m — d 

4. NON-LINEAR OPERATORS 

Our objective in this section is to show that the results of the two preceding 
sections can be applied to prove a structure theorem for a class of non-linear differ- 
ential operators admitting invariant polynomial subspaces. We shall see that these 
results complement the structure theorems for operators preserving simplicial mod- 
ules which were proved in [14]. 

In dealing with non-linear operators it is convenient to identify differential op- 
erators with functions on jet space. To that end, let 

J''(M) = M X W+^ 

denote the bundle of r-jets of smooth maps from R to M. The r-th prolongation of 
a smooth, real- valued function /(x) is a section of J"^ , namely 

pTXf) = if,Df,D^f,...,DJ) 

Thus, the action of an operator on a function of x is the same thing as the compo- 
sition of a function of the jet variables with the prolongation: 

T[/]-ropr,(/). 

We introduce the standard jet coordinates x, uq = u,ui = Ux,U2 = Uxx, ■ ■ ■ ,Ur on 
J^ so that 

D^[f]=u,opr,if)- 
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Henceforth, we fix n. By Propositions |H1 and El all linear operators of order r 
are expressed uniquely as polynomial linear combinations of the maximal deficiency 
operators Krj- Thus, a linear operator Krj maps Vn to Vn-r] a quadratically non- 
linear KriKrj maps Vn to 'p2(n-r)j a cubically nonlinear KriKrjKrk rnaps Vn to 
Vz{n-r)i etc. This implies: 

Proposition 11. Every operator ( linear or non-linear ) of order r can be uniquely 
expressed as 

(25) T := p{x)+^^Pi{x)Kri+^^Ptj{x)KriKrj+ ^ Pijk{x)KriKrjKrk^ , 

i i<j i<j<k 

The operator in question will map Vn to Vn */ o,nd only if 

(26) P eVn, Pi &Vr, Pij e V2r-n, Pi]k G V3r~2n, ■ ■ ■ 

The above proposition has the following obvious consequence: 

Corollary 1. An operator of order r < ^, mapping Vn to Vn is necessarily linear. 
If r < ^, then the operator will have at most quadratic non-linearities. 

Conversely, the above proposition can be used to bound the degree of any poly- 
nomial space which can be left invariant by a non-linear operator. We have, for 
example: 

Corollary 2. A second-order operator will preserve a polynomial space of degree 
at most four. 

The extension of these results to multivariate differential operators acting on 
simplicial polynomial modules is straightforward. We are interested in writing 
down all non-linear differential operators of order r that preserve 

(27) Vn = {x\' ...x'JJ \ii + --- + iN< n}. 

We recall that for fixed order r the maximum deficiency that an operator can attain 
is r, which is achieved by any of the following maximum deficiency operators: 

(28) Kj e Cr,r, =^ Kj - (n - \J\ - E),_ij\ D\ 

(29) J = {ji, ...,JN\ji + ---+jN < r}. 

The extension of Proposition^Jto the case of multivariate differential operators 
and polynomial modules is straightforward by just substituting the simple indices 
i,j into multi-indices /, J, each of which can assume dim^^.r = ( ^^) different 
values. Similarly, the two Corollaries hold verbatim in the multivariate case. 
Example 2. Write down all second order operators with quadratic non-linearities 
that map V4 into itself. A direct application of Proposition II II allows to write: 

(30) T[u]= Y. PvK2iK2, 

0<i<j<2 

where Pij — pji are constants, and where 

K20 = 6uo — 3xui + i x^U2, 

K21 = 3mi - XU2, 
K22 = U2 
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The resulting quadratic combination is the following six parameter family of oper- 
ators: 

T = poou\+Pol{-Xu\+iuiU2) +PQ2{x'^u\~&XUlU2 + l2u2Uo) 
+ Pll {x U2 + QUi — QXU1U2) 

+ P12 {—x^u^ + 9x^uiU2 — 6x{2u2Uo + 3ui) + 36uiUo)) 

+ P22 {j x*U2 — 3x^uiU2 + 3x^{2u2Uo + Suf) — 36xuiUo + 36ug) 

If we are interested in finding only autonomous non-linear equations we need to 
consider the translation-invariant subfamily of the above operators, i.e. operators 
where the variable x does not appear in the coefficients. Imposing this condition 
leads to 

Pio = P12 = P22 = 0; Pii+P2a = 0, 
so that the only autonomous second order non-linear operators which preserve the 
space V4 are: 

(31) {ul^ , AuxxU - 3ul}. 

The first of these two operators is easily seen to map V4 into itself. The second 
operator is more interesting. Each of the non-linear terms u^ and UxxU transform 
a 4"^-degree polynomial into a 6*'^-degree polynomial. However, the linear combi- 
nation 4uxxU — 3u^ cancels the coefficients of degree 6 and 5, and hence defines a 
map from V4 to itself (see also [15] and [19]). 

We investigate further in the analysis of autonomous non-linear operators with 
invariant polynomial subspaces in the following section. 

5. The algebra of polynomially non-linear operators 

In this section we continue the study non-linear operators that are polynomial 
in the function and its first n derivatives. First, let us fix some notation. Let 

T = M[.T,Uo,. .. ,Uri] 

denote the commutative algebra of non-linear operators that can be expressed as 
polynomials over M in a; and the derivatives Uj . Multiplication in this algebra is by 
pointwise multiplication, rather than operator composition. We grade this algebra 
by total degree in the Uj variables: 

00 

/=0 

where 

Te = spanja;-' Ui^u^^ ■■■Ui, | < j < 00}. 

We will refer to the integer d as an operator's degree of non-linearity. Thus, 7^ 
is the vector space of linear operators, 7^ the vector space of quadratically non- 
linear operators, etc. The vector space Tq is the space of constant operators. For 
example, the operator x maps all of Vn to x. Thus, Tq is a subalgebra of T, which 
is isomorphic to the polynomial algebra in the variable x. All the other 7^, i? > 1 
are merely subspaces of T, not subalgebras. 

We further grade each 7^ according to the following monomial weighting scheme: 

(32) wt(Mi) =11 — 1, wt(a;) = 1. 
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Thus, 

OO 

(33) T, = 0r,,fe, 

k=Q 

where 

e 

(34) Ti^k = spanix-' Ui^Ut^ ■■ -Ui^ | j + n^ - ^ is = fc} 

s=l 

is the subspace generated by monomials having weight k. We will refer to the 
integer n — k as the monomial deficiency. If M = x^Ui^Ui^ ■ ■ -Ui^ is an operator 
monomial with 

i 

(35) fc = wt(M) = j + n£-^is, 

s=l 

then, in accord with the above-introduced meaning of deficiency, M maps Vn into 
■Pfe, but not into Vk-i- 

In other words, every operator T E T admits the unique decomposition 

(36) r-^r,,fc, 

e,k 
where 

k=j+ne—(ii+i2-\ \-it) 

and where the sum is taken over finitely many values of £ and k. For generic values 
of the coefhcients Cji-^,,,i^, the operator Ti^k has deficiency n — k. However, for 
certain very specific values of the coefficients, the actual deficiency is greater than 
the monomial deficiency. This is so because in the linear combination H37|l there 
might occur some cancellations in the terms of highest degree. We see for instance 
in Example 2 that relative to 7^4, the operator 

4:UxxU — 3u^ = 4:U2Uq — 3u1 

has monomial deficiency 4 — fc = 4 — 8 + 2 = —2. However, the actual deficiency of 
this operator is zero. 

Next, we describe generators for T that will allow us to precisely determine the 
deficiency of an operator. Following Proposition EJ let us re-introduce the linear, 
n^^-order operators of maximal deficiency: 

n-j i 

(38) Vj = Knj = ^{-'^y -j- Ui+,j j = 0,...,n. 

i=0 
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Thus, 






(39) 


Vn = W„, 




(40) 


Wn-1 = Wn-1 - 


-XUn 



Wn-2 = Un-2 - XUn-1 + \ X^Un-2, 

Wn-3 = ""n-S - XUn-2 + 5 a;^U„_3 - \ X^W„_4, 



These operators transform all elements of Vn into a constant. In other words, the 
Vj are the operator duals [14] to the monomial basis of Vn'- 

Vj[x''/k^] = 5j^-, fc = 0,l,...,n. 

More, generally let us write 

(41) Uj [t) = ^ w,+j - j = 0, . . . , n. 

In this way, 

V, =uA-x). 



Proposition 12. We have 

(42) w,(s + i) = ^u,+j(s)^. 



E 

i=0 



Proof We have 



u 






'{t) ^Uj+i{t), j = 0, l,...,n- 1, 

Un'it) = 0, 

Hence, Uj i— > Uj{t) defines a 1-parameter transformation group of T. The desired 
result follows immediately. D 

Now, we can invert the relations (|38|) , and express the Uj in terms of the Vj . 

Proposition 13. For j = 0, . . . , n, we have 

n-j j 

(43) Wj=^w,+j— , 



i=0 



Proof We apply 142(1 with s = —x and t = x. D 

Proposition 1131 shows that a;, wq, . . . , fn freely generate the algebra T. The rela- 
tions (|38|l and (|43|) are homogeneous relative to the weights H32II . Hence, setting 

(44) vft{vj)^n-j, 

we recover the grading by monomial deficiency relative to this basis. We now deepen 
the grading by defining 

i 
Te,k,m = spanja;" Vi^Vi^^ • ■ -Vi^ \ m + n£ -\Js = k}, 

s=l 
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SO that 



%,k — ^j3 ■^^'=^' 



Proposition 14. The elements ofTi^k,m have monomial deficiency n — k and actual 
deficiency n — m. Consequently, every operator T € R[x, uq, . . . , Un] has deficiency 
n — m, where m is the x-degree of the polynomial 

T = Q(a;,t)o,. ■ ■ ,Vn) 
that expresses T relative to the x, Vj basis. 

Proof Since the relations (|38|) (|43|) are homogeneous relative to the weighting 
scheme (|32ll (|44() . the elements of T£^k,m have weight k, and hence have monomial 
deficiency n~k. Since the operators Vj map all polynomial to constants, an operator 
Q{x,vo, . . . ,f„) e T, having m as its x-degree, maps Vn to Vm, but not to Vm-i- 
Therefore such an operator has deficiency n — m. D 

The key application of the above grading has to do with the decomposition 
of an operator according to monomial deficiency. Our analysis would be greatly 
simplified if we could be certain that the decomposition of a non-linear operator 
T according to monomial deficiency respects the actual deficiency. In other words, 
when considering operators of a fixed deficiency, no generality is lost by considering 
operators that are homogeneous in degree of non-linearity and monomial deficiency. 

Corollary 3. Let T be a non-linear operator whose deficiency is n ^ m or more, 
i.e., T maps Vn into Vm- Let Tg^k be the summands of the decomposition of T 
according to degree of non-linearity and monomial defiency as per (|36f) (|37|l . Then, 
each Ti^k also maps Vn into Vm- 

6. Autonomous, non-linear operators 
Our main focus in this section is the subalgebra 

(45) ^ = M[mo,--- ,u«]cr 

of translation-invariant non-linear opearators. The subalgebra inherits the bi- 
grading relative to degree of non-linearity and monomial deficiency, with 

ctt nt 

(46) ^ = 00^,^,, 

t=0 k=0 

where 

(47) A,fc = .4 n 7^,fc = span{ui^Uj2 • • • Uj^ | n£ - ^ z^ = k} 

s=l 

Our key result in this section is the characterization of the deficiency of autonomous 
operators. In other words, we describe ^fl Ti^k,m- 

The obvious approach to construct non-linear operators of deficiency m would 
be to write a generic polynomial p(x) S Vn with indeterminate coefficients, act on 
it by l|^ where the degree of the possible non-linearities is bounded by Proposition 
nil and impose that the coefficients of all the terms in x^ for j > n — m vanish. 
However, based on the useful concept of deficiency we choose here to adopt a 
somewhat different approach. 
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We seek operators that are both translation-invariant and that have maximum 
deficiency. To this end, we define 

(48) C = ^^^; 

(49) /„_,={tj(-e)-^(-l)'"z+,|, j = 0,...,n, 

i=0 

and note that 

(50) /o = w„, /i = 0. 

The jet space function ^ is only defined on open neighborhood ?i„ 7^ of JT". 
Thus, the operators Ij are defined only for elements of V^ — VnX Vn-i, the set 
of polynomials of degree exactly equal to n. We can still speak of the deficiency of 
such operator, but this has to be understand in terms of V^ rather than Vn- 

Proposition 15. The non-linear operators I2, I3 ■ ■ ■ , In are translation-invariant, 
and have maximum deficiency n. In other words, these autonomous, nonlinear 
operators transform every n*"^ degree polynomial into a constant. 

Proof By definition, 

«J = Uj{~x), In-j = Uj{-£,). 

As well, 

(51) ^ = ^ + !!iizl. 

Hence, Proposition 1121 implies that in addition to H49() we also have 

(52) ^"-.-E(-i)'^-^+.R-" 

Hence, the operators Ij are polynomials of the Vj divided by a certain power of 
Vn = Un- Therefore, these operators are both translation invariant and of maximal 
deficiency. D 

We can also invert the relations H49|l . and express the Uj in terms of the Ij. 



Proposition 16. For j - 


= 0, . . . , n, we have 


(53) 


2 — 


('54^ 





i=0 



[n- i -i)V 



Proof We apply ^^ with s = -S, and i = ^. D 

Thus, relations (|I^ ifK!^ tell us that the operators /q, ^jh ■ ■ ■ ,In also generate the 
algebra of autonomous operators. These relations are homogeneous with respect to 
the monomial weights defined in l|S^ . Hence, setting 

(55) wt(/,) = J, wt(0 = 1 

we recover the grading by monomial deficiency relative to this basis. 



DIFFERENTIAL OPERATORS ADMITTING INVARIANT POLYNOMIAL SUBSPACES 15 

Unfortunately, the operators Ij , ^ are not polynomials in the Uj , and hence are 
not elements of A. Let us therefore consider the larger algebra 

(56) i = R[/o,e,/2,/3, •••,/„] 

generated by ^ and the autonomous operators of maximal deficiency. Thanks to 
(|5;-i|l we know that A C A, but this inclusion is strict. 
We now deepen the grading of A by defining 

Ai^k,m = span{^™ hi- ■ -lit I m + ii H Vii = k}, 

so that 



oo nl 



-4 = 000^ 



e=o fc=0 m=0 

Proposition 17. We have Ai,k.m = -4 fl ?£ j, ,„, i.e., the elements of Ae^k,m have 
monomial deficiency n ~ k and actual deficiency n — m. 

Proof This follows from l|51|) (|52|) , and the fact that the operators Vj have maximal 
deficiency. D 

In other words, the deficiency of an autonomous operator is n minus the ^-degree of 
the polynomial that expresses that operator relative to the ^, Ij basis. We are left 
with the question of the nature of the inclusion of A in A. In other words, which 
polynomials in ^, Ij define true polynomial operators. 

Theorem 1. Let T — P{x, uo, . . . , Un) £ T be a non-linear operator, and let T = 
Q{x,vo, . . . ,Vn) be the expression of this operator relative to the non- autonomous 
generators x,Vj. Then T is autonomous, i.e. P is independent of the variable x if 
and only if 

T = Q(e, /„,..., /2,0,/o). 
In this case, the deficiency ofT is equal to n minus the £,-degree of the polynomial 

Proof By Propositions II 31 and II 61 we have 

i=0 ' i=0 

where, as we noted before, /i = 0, and /q = w„ = u„. D 

Example 3. Let us recast the analysis began in Example 2 in terms of the above 
operator bases. Relations (|43|) and H53|) take the form 

U4 = W4 = /o ; 

U3^V3 + XV4 ^CIo, 

U2 = V2+ XV3 + i a;^U4 = h + ^ S."^ h 

ui =vi+ XV2 + I x'^V3 + i x^Vi = h + S,l2 + ^ £.^Io; 

Uo = Wo + xvi + \ x^V2 + i x^v-i + jg x'^Vi =74 + ^/3 + 1 ^^/2 + jg Sf'h- 

Our goal is to write down all autonomous operators with quadratic non-linearity 
and zero deficiency, i.e., operators that map Va into itself. By Proposition 1171 
and the above relations, we are obliged to consider polynomials that are quadratic 
in lo, Ii, l2, I4 and that have degree 4 or less in the ^ variable. The question 
is: which operators of such form are quadratic in mo,mi, 1*2,^3,114? By Corollary 
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no generality is lost by considering operators of a fixed monomial deficiency. 
Evidently, if the monomial deficiency is or more, the operator will preserve Vi- 
Hence, we must consider operators having monomial deficiency —4, —3,— 2,-1, 
which corresponds to fc = 8, 7, 6, 5, respectively. The most general operator having 
A: = 8 is a multiple of the monomial (uo)^- Such an operator will have actual 
deficiency of —4 as well. Similarly reasoning holds for operators with fc = 7; these 
are multiples of uqUi. Let us consider the case fc = 6. The ansatz is now 

Co2 U0U2 + Cii uj ^ ^44^^*^°^ '^ ^'^"^ ''^ a^^ + ^(^^02 + 4Cii) v^v^ x^ + l.o.t. 

= Yii(3Co2 + 4Cn) lie + Oe' + l.o.t. 



Hence, such an operator preserves Va if and only if is a multiple of AuqU2 — 3wf , as 
has already been proven by other methods. 

Finally, let us consider the case fc = 5. Now, the ansatz is 

1 5 

Co3 "0^3 + Ci2 UiU2 = ^(C03 + 2C12) vl X^ + ^(Co3 + 2C12) VsV^X^ + l.o.t. 

= ^ (Co3 + 2C12) II e + Oe + l.o.t. 

Hence, such an operator preserves Vi if and only if is a multiple of 2uoU3 — uiM2. In- 
deed, because /i = 0, the above calculation shows that this operator has deficiency 
f , i.e. it maps Vi into V3. 

Autonomous operators with quadratic non-linearities. We restrict from 
here on to 7^, the vector space of operators with homogeneously quadratic non- 
linearity. We present a complete characterization of such operators in terms of de- 
ficiency, thereby extending the results of Svirshchevskii [19] and Galaktionov [15]. 
Our analysis can be extended to operators with higher non-linearities, but this shall 
be treated elsewhere. 
Following |j22|), we let 

n 

Q= ^ Qk 

be the linear space of quadratic autnomous non-linear operators up to order n, 
graded according to monomial deficiency. Thus, 

(57) Q = A2 = spa,n{uiUj | < i < j < n}, 

^n-l-2^ 



(58) dimQ 2 

(59) Qk = A2,n-k = spa.n{uiUj \ i + j = n + k, <i,j <n}, 



(60) dim Qk 



1, 



where and [-J denotes the floor function. 

Each of the elements of the above basis of Qk has a different order r, the minimum 
and maximum orders for each fc being: 

(61) r^^ - [^1 , rW, = min(n, n + fc) 
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All of the above basis elements have the same deficiency k, i.e. they map Vn into 
Vn-k- However, Example 2 and Example 3 show that within each Qk there can be 
elements whose deficiency is greater than k. 

Remarkably, each Qk possesses an adapted basis that further grades it according 
to operator order and deficiency. In this regard, for r^l^ < ?" < ?'inax, let us 
introduce the operators 

(62) Qk,r= y](-l)M* , A\ ]U,Uri+k-i, 

^ — ' \n — r+l/Vn — r/ 

i—k-\-n—r 

We are now ready to state the main result of this Section: 

Theorem 2. The operator Qk.r has monomial deficiency k, order r, and deficiency 

(63) m(fc,r) = fc + 2(r-r|„l). 

Furthermore, {Qk,r \ TmL < ^ < rinlx} forms a basis of Qk- 

This result is similar to that proved in Proposition 0] for linear operators: up 
to a scalar multiple there is only one quadratic autonomous operator with a given 
order r, monomial deficiency k, and deficiency m(k,r). In Appendix A we show 
the explicit form of these operators for n = 4, 5 and 6. 

We will prove the theorem by finding a generating function for the operators 
Qk,r- This requires some new notation. For a bivariate formal series, 

p{z, w) = ^ p^j z'w^ , 
i,j>0 

let us define 

(64) g{p{z,w)} = Y,P^.^z' 

i>0 

to be the series formed from terms where the two variables have equal exponents. 
We will also adopt the convention that 

/n\ f q(q-1)---(q-^ + 1) . . r, 

(65) (")= il ' '^' 

VV [o i<0 

Proof [ Proof of Theorem|21] We extend the definition (|62|l of Qk.r to all < r < n 
by setting 

min(«^+fc) ./i^k+l\/n-i\ 

(66) Qk.r= y^ (-l)M , -, ]uiUn+k-i- 

i=max(0,fc) ^ / \ / 



18 DAVID GOMEZ-ULLATE, NIKY KAMRAN, AND ROBERT MILSON 

We can now form a generating function for Qk,r as follows: 

(67) 

n n 

^n — r _i_n — k 



Q{z,t)^ E E^^-^""'*" 



k——n r— 

i,j=Or = \ / \ 



n n / \ s 

(68) =EE(-i)1 ]z''g{z{i+zr-\i+wr-'+'}i,i,^t^"-^-', 

where s — 2n — i ~ j — p — q, and where we used relation (|54|) to replace the Uj 
with the Ij , the autonomous generators of maximal deficiency defined by (|49|) . Let 
us extend the definition of Ij to all j > by setting 



(69) /, = E(-l)^"n+^-A, 



and agreeing that Uj — for j < 0. However, we must note that the newly 
defined autonomous operators Ij , j > n are no longer maximal deficiency operators. 
Interchanging the summation order in (|68|l and re-indexing with 

i —>■ n — p — i, j —>^ n — q ^ j 

we obtain 

n n-pn-q / \ 

Q{z,t)^ E EE(-i)'L_',-h"'^{^(i+^)""'(i+H""'^'>^''^'^^ *'""'"' 

p,g=0 i=0 j=0 \ / 

n n-pn-q / \ 

= E EE(-i)M •^"'^{^(l + ^)'^"(l + ™)'^''''>^''^'^^*'^''^^ 

with s — i + j henceforth. Relation (|54|l holds for j < 0, thanks to the extended 
definition H()9|) oi Ij, j > n. Hence, 

2n 2n — p — q 2n — p — q — i / \ ^ 

Q(^,i)=E E E (-i)l'U"^G{^(i + ^)''+'(i + ^)^+^+'}v,^i''+'+^ 

p,g=0 i=0 j = \ / 

2n 2n-p-q " / \ ;^s 

= E E E(-i)""''"M Jp~'^{^(i + ^)''^"(i + '^)'^'^'""}^j'^^7[*''^'^' 



(™) = E E {-^T~^z'^Q{4i + zf{i + wY^\w-zy}ipi,^t^+''+' 

p,q^Q s— 

For every positive integer s, let us introduce the generating function 

oo 

(71) 0(p, q, s; z) = g{{l + z)P(l + ivfiiD - zY) - E -^P^s^P.- ^''' 

where 

'^^' *— i:<-"'i')(2.-V,-)C-. 
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Thus, 

(73) 4'p,q,p,s =0, S > 2/9, 
(Pp,q,p,2p-1 = iq-p)i 

s-©)rr)--c;:r 

Before continuing, let us note the foUowing properties of this function: 

(74) (j){p, q, s; z) = {-iy(j){q,p, s; z) 

(75) 0(p, g, s + 1, 2) = (l){p, q + l,s;z)- (l){p + 1, 9, s; z) 
This function is relevant to our proof, because 

(76) g{z{l + zfil + w)«+i {w - zf} ^ 0(p + 1, g + 1, s; z) - 0(p, g + 1, s; z). 
Hence, 

n n 2n 

(77) Q(^,i)= E EE^fc.' 






;,r.s , -2^ ^ 



k——n r—0 s—0 

where 

(78) Qk,r,s = Qk,n~p,s = / (~1) (0p+l,g + l,p+l,s ^ 0p,g+l,p+l,s) -^p-^q ; 

p+q^Ti — fc — s 
n 

(79) ^ Qfe..,.^"-'- = ^ {-!)"-" z-\cb{p+l,q+l,s;z)-^{p,q+l,s;z))lpl,. 

r— p+q^n— fc— s 

Let m(k,r) denote the deficiency of Qk,r- By Proposition^ m{k,r) is equal to n 
minus the largest value of s for which Qk,r,s 7^ 0. By l|73(l and (|78|l . we know that 
Qk,n~p,s = for s > 2p + 2. For s = 2/9 + 2, we have 

^ w, -r nP+i/^/^2p + 2\ /2p + 2\\ 

<?>p+l,g+l,p+l,2p+2 - ?'p,9+l,p+l,2p+2 - (-tj M i]^l^l _|_l))~' 

and hence, Qk,n-p.s — for s = 2p + 2. 

The analysis now breaks up into two cases. Suppose that n — fc is even. For 
s = 2p + 1 we again use ^TM and (|78|l to obtain 

gfc,„-p,2p+i = E (-1)"-^+'^+^ {^^ ^ ^) ipiq - 0, 

p+g=n-fc-2p-l \ P J 

because of p, q symmetry. For s = 2p, we have 

Qfe,„-P.2p- 1^ (-1) 2(p+l) (,JV,^o. 

p+g=n-fe-2p ^'^ ^ ^ '^ ^ 

Therefore, if ?i — fc is even, we have m{k, r) = n — 2p = 2r — n. This agrees with 
Next, suppose that n — fc is odd. Now for s = 2p + 1 we have 

Ofc,„-p,2p+i = E (-1)"-^+"+' (^^ ^ ^) ipiq ^ 0, 

p+g=n-fc-2p-l \ P / 
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because now p and q have the same parity. Therefore, iin — k is odd, then m(fc, r) — 
71 — 2p — l = 2r — n — 1. Again, this agrees with 1)63(1 . D 

7. Separation of variables in non-linear evolution equations 

One of the main appUcations of the resuhs of this paper Hes in the method of 
separation of variables in non-Unear evolution equations. Let u £ Vn and T be a 
non-linear differential operator with deficiency tti > 0. This means that TVn C Vn, 
or more explicitly 



^ 



/ ^PiX 



^^MPl,---,Pn)x'' 



(80) T 

=0 J i=0 

The non-linear evolution equation 

(81) Ut^T[u]^p{x,U,Ui,...,Un) 

admits separable solutions of the form 

n 

(82) u{x,t)=Y,V^{t)x\ 

4 = 

where the functions (p{t) satisfy the following system of first order ordinary differ- 
ential equations: 

(83) (pi = fi{ipi,...,(pn), i = 0, ...,n. 

These results extend immediately to the multivariate case. This technique has been 
used by King to find new exact multidimensional solutions of non-linear diffusion 
equations [16, 18] and exact solutions of high order thin film equations [17]. In a 
more mathematical context, Galaktionov [15] and Svirshchevskii [19] have analyzed 
non-linear operators that preserve low dimensional spaces spanned by polynomials 
and trigonometric functions. 
Example 4 Consider the following non-linear evolution equation for u = u{t, x): 

r.( 5 45 2 

(84) Ut^" \U U^^r^x - -jUx U^xx + ^ "x 

This equation has solutions in Vg, of the form 

8 

(85) M(x,t)=^(^,(i)a;\ 

i=Q 

where the functions ipi{t) satisfy the following first order system: 

ipa = 45<y92 ~ 105<y9i(p3 + 168 (/3o'^4 

ipi = 60 ip2(p3 - 252 tpitpi + 840 ipo(p5 

(p2 — 90tpl — 132Lp2'P4 — 210ifiip5 + 2520Lpoipe 

ips — 108 (ps(p4 — 360 ip2f5 + 4:20 (fiifa + 5880 ^0^7 

ip4 = 108(^4 — 135 1^931^95 — 330 (y92<y26 + 2205 iy9iiy97 -I- 11760 (y9o¥'8 

ip5 — 108 (p4ip5 — 360 (fi^ipQ + A20 (p2'P7 + 5880 (fiiifs 

ipQ = 90(y9g — 132 (ya4(y96 — 210(/J3(^7 -I- 2520(/32'y38 

ipr — 60 ip^ipe — 252 ip4(pj + 840 tpsips 

(fs — 45 (y9g — 105 (p5ip7 + 168 ip4<fs 
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Appendix A 

In this Appendix we list the speciahzed autonomous operator basis Qk.r de- 
scribed in Theorem 1. Recall that each Qk,r has monomial deficiency fc, order 
r, and deficiency m = m{k,r). The following tables for n = 4,5,6 display the 
non-linear operators with {fc,r, to} in the ranges 

. + k' 



< k < 



n, 



< r < min(n, n + k), m ^ k + 2(r — rmin) 



Recall that an operator with deficiency to- maps Vn to Vn-rm but not Vn-m-i- 
Hence all operators with to, > map Vn to Vn and can be used to construct an 
evolution equation solvable by non-linear separation of variables. 

Table 1 . Quadratic autonomous operators acting on P4 



k 






Qk,r 








to 






-4 


2 












-4 






-3 


UqUi 












-3 






-2 


ul 


3ul- 


- 4:U0U2 








-2 







-1 


UIU2 


U1U2 - 


- '^UqUs 








-1 


1 







ul 


2u\- 


- 3miU3 


ul- 


- 2miU3 


+ 2uoU4 





2 


4 


1 


U2U3 


U2U3 - 


- 3uiU4 








1 


3 




2 


ul 


uh 


- 2U2U4 








2 


4 




3 


U3U4 












3 






4 


ul 












4 







Table 2. Quadratic autonomous operators acting on V5 



k 






Qk,r 




m 






-5 


2 

"0 








-5 






-4 


UqUi 








-4 






-3 


ul 


Au\- 


- 5uoU2 




-3 


-1 




-2 


UIU2 


'iU\U2 - 


- 5uoU3 




-2 







-1 


ul 


3ul - 


- 'iuiU3 


9ul — I6U1U3 + 10^0^4 


-1 


1 


3 





U2U3 


U2U3 - 


- 2uiU4 


U2U3 — 3uiU4 + 5uoU5 





2 


4 


1 


ul 


2ul- 


- 3U2U4 


ul — 2U2U4 + 2UlUr, 


1 


3 


5 


2 


U3U4 


U3U4 - 


- SU2U5 




2 


4 




3 


ul 


ul - 


- 2U3U5 




3 


5 




4 


U4U5 








4 






5 


ul 








5 
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Table 3. Quadratic autonomous operators acting on Pq 



k 






Qk,r 










m 






-6 


ul 














-6 






-5 


UO Ul 














-5 






-4 


ul 


5ui - 


- 6uoU2 










-4 


-2 




-3 


U\U2 


2U\U2 - 


- 3uoU3 










-3 


-1 




-2 


ul 


-iul - 


- 5uius 


6U2 — lO-UlUs + 5uoU4 








-2 





2 


-1 


U2U3 


3U2M3 - 


- 5uiU4 


2U2U3 — 5uiU4 + 5uoUs 








-1 


1 


3 





ui 


3«i- 


- 4U2U4 


9ul — I6U2U4 + lOuius ul - 


~ 2U2U4 


+ 2uiU5 - 


- 2uoU6 





2 


4 6 


1 


U3U4, 


U3U4, - 


- 2U2U5 


U3U4 — 3W2M5 + 5uiU6 








1 


3 


5 


2 


ul 


2ul- 


- SU3U5 


u\ — 2M3M5 + 2U2U(, 








2 


4 


6 


3 


U4U5 


U4U5 - 


- 3U3U6 










3 


5 




4 


ul 


ul' 


- 2U4U% 










4 


6 




5 


U5Ue 














5 






6 


ui 














6 
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